The ideal relativistic spinning gas: polarization and spectra 



F. Becattini 

Universita di Firenze and INFN Sezione di Firenze, Florence, Italy 



i — \ 



F. Piccinini 

INFN Sezione di Pavia, Pavia, 

We study the physics of the ideal relativistic rotating gas at thermodynamical equilibrium and 
provide analytical expressions of the momentum spectra and polarization vector for the case of 
massive particles with spin 1/2 and 1. We show that the finite angular momentum J entails an 
anisotropy in momentum spectra, with particles emitted orthogonally to J having, on average, a 
larger momentum than along its direction. Unlike in the non-relativistic case, the proper polarization 
vector turns out not to be aligned with the total angular momentum with a non-trivial momentum 
dependence. 

I. INTRODUCTION 

The thermodynamics of rotating systems is a subject of interest mainly in astrophysics, where spinning objects are 
compact stars. Also in nuclear physics, there have been noteworthy applications to the problem of multifragmentation 
[HQ. In most those applications, either the non-relativistic approximation is used or the spin of particles is neglected, 
so that a full description of the relativistic case including spin degrees of freedom is still missing. However, recently, 
the microcanonical and grand-canonical partition function of an ideal relativistic quantum gas of particles with spin 
have been calculated Q enforcing angular momentum conservation. Taking advantage of these results, in this paper 
we work out analytical expressions of the spectra and polarization of particles in a relativistic rotating gas with large 
(in h units) angular momentum. From a phenomenological point of view, these calculations might be of interest for the 
physics of relativistic heavy ion collisions, where the formation of a system with a large intrinsic angular momentum 
has been envisaged 0] . We will confine ourselves to the case of Boltzmann statistics, leaving the quantum statistics 
case to future work. 

The paper is organized as follows: in Sect. [IT] we discuss the statistical mechanics of an ideal relativistic gas with 
fixed, large angular momentum and show the equivalence with a rigidly rotating system; in Sect. IIIII we analyze 
more in detail the relation between angular momentum and angular velocity in the limit of low rotational speed; in 
Sect. HVl we calculate the inclusive momentum spectra of particles; finally, in Sect. |V] we obtain the expressions of the 
polarization vectors of massive particles with spin 1/2 and 1. 

Notation 

In this paper we adopt the natural units, with h — c = K = 1. Space-time linear transformations (translations, 
rotations, boosts) and SL(2,C) transformations are written in serif font, e.g. R, L. Operators in Hilbert space will be 
denoted by an upper hat, e.g. R. Unit vectors will be denoted with a smaller hat e.g. p. Even though the notation 
is unambiguous, we will make explicit mention of either possibility whenever confusion may arise. 



II. ROTATING RELATIVISTIC GAS AT EQUILIBRIUM 



In general, by rotating thermodynamical system we mean a system with fixed angular momentum in its rest frame. 
This is a good definition also in the relativistic case. It is well known that a classical system with non-vanishing 
intrinsic angular momentum can be at thermodynamical equilibrium only if the rotation is rigid ||, that is if a 
constant vector uj exists such that the local collective velocity v is: 

V = U) X X (1) 

We will show that this conclusion holds for the relativistic gas, with the obvious inequality: 

\\uj x xj| < 1 (2) 

In this section, we will develop a proof based on an ab initio statistical mechanics calculation while a proof based on 
an extension of Landau's argument for relativistic systems is given in Appendix B. 
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In ref. [|| , one of us derived the full expression of the microcanonical partition function, as well as its grand-canonical 
limit for large volumes, of a multi-species ideal relativistic gas with fixed angular momentum. In relativistic quantum 
mechanics, fixing the angular momentum means projecting onto irreducible states of the orthochronous Poincare 
group 10(1,3)^ which, in fact, is a projection onto angular momentum states in the rest frame, where the total linear 
momentum vanishes. If, moreover, the intrinsic angular momentum is large (in Ti units), it can be treated as a classical 
vector J and the grand-canonical partition function reads 0] : 



Zj = tr{cxp[(-tf + M Q)/T]PjPv4 
2J+ 1 



d"0 <> 

</>|<ir 



i<p-J 



exp 



J d 3 x J d 3 pe-^/ T e-^ ( xx P)tr^^(R < ^(< ? !>)) 



(3) 



where T is the temperature and V the volume; £j is the energy, Xj the fugacity and Sj the spin of the jth particle 
species; D s ^ (R^((f>)) is the matrix of the irreducible representation Sj of SU(2) transformation R with axis and 

angle 4> — \\4>\\. The operator H is the hamiltonian, Q is a generic conserved charge operator and \i the relevant 
chemical potential; the operator Pj is the projector onto a fixed angular momentum while: 



p v = ^2\h v ){h l 



is the projector onto localized states \hv), which form a complete set of quantum states for the system in the finite 
region V. This projector is needed in eq. © for the trace operation to be a properly defined one, i.e. involving a 
basis of the full Hilbert space [|| . 

The formula (|3|) applies to a system with large V and J and the integrand function turns out to be significant only 
in a region where <f> — \\4>\\ <C 1. Since both J and V are large, it is then possible to make a saddle-point expansion of 
the integral ([3]). To do this, we first have to continue the integration variables to the complex plane and, in view of 
the spherical symmetry of the domain, the obvious choice is to take the spherical coordinates of the vector <f>, namely 
its magnitude (j> along with its polar and azimuthal angles. Then, we have to solve the complex vector equation: 



V fi 



* ' J + E 7^3 / d3x / d3 P ^ ( R <^)) 



e £j/r e -i<Mxxp) 



(4) 



leading to: 



J = // 3x / d 3 P tr^(R^(0))(xxp)e-^/ T e-^( xx P) 



(27T) 

J 

+ E (^// 3x / d3 P [^trD^R^)) 
= L(0) + S(^») 



(5) 



the reason for naming the two integral terms as L(</>) and S(<p) will become clear later on. We first observe that, if 
the equation (gj) is solved by cp, <f>, then also —<p*,(p is a solution. This can be easily checked by taking into account 
that the trace depends only on <fi and can be written as a sum of exponentials exp[m^>], being n an integer. Therefore, 
we will look for one solution enforcing (f> — — (f>* and <j> = <p , that is with a real unit vector cj) and an imaginary 
magnitude <j). 

Assuming that one solution exists and defining the real vector uj = —iT<fi, we can write the grand-canonical partition 
function at the lowest order of the saddle-point expansion as: 



Zj cx exp[— us ■ J/T] exp 



J d 3 x |d 3 ptr^(R^(zc i ;/r))e-^/ T e-( xx P)/ T 



(6) 



The function: 



exp 



y -\ [ d 3 x I d 3 ptr J D^(R i (^/r))e-^/ T e w ( xxp )/ T 
. (2ti") 3 Jv J 



(7) 
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is the partition function of an ideal relativistic gas (in the Boltzmann limit) rotating with an angular velocity u>. This 
can be shown explicitely for the spinless case by dividing the gas into small cells with volume A 3 x and uniform velocity 
ffl and calculating the relevant grand-canonical partition function in the formalism of relativistic thermodynamics 



C ccU = exp 



E 



(2tt) 3 



A 3 x / d 3 p e 



-P-p 



where (3 is the temperature four-vector: 

(3 = 



= (l,v) = 



1 



:(1,W X x) 



(8) 



(9) 



To being the local temperature, i.e. the temperature measured by a thermometer moving along with the cell. Note 
that with this equation, we are tacitly assuming that the partition function, which is a Lorentz-invariant quantity, in 
the accelerated rotating cell is the same as that in an inertial frame with its instantaneous origin and velocity 1 . This 
is consistent with the locality hypothesis 0] and the transformation law of energy in a rotating frame 0: 



E' = j(E - u> • L) = i(E - u> ■ (x X P)) 



being 7 = l/yl — ffw x x|| 2 . Therefore, with the identification: 



t = Wi-« 2 = Wi - ¥> x x l! 2 

the partition function of the cell becomes: 



Z ce n = exp 



^(^ A3x / d3pe "" /TeP( " XX)/T 



(10) 



(11) 



(12) 



The total partition function of the supposedly independent cells can now be obtained by multiplying expressions JT 
for all cells and going to the limit of infinitesimal cells, implying: 



exp 



E 



A, 



'2tt) 3 



d 3 x / d 3 p e 



-e 3 /r e ^-(xxp)/T 



(13) 



which is the ([7|) for spinless particles. 

Therefore, unlike in the non-relativistic case, the proper temperature To m a spinning relativistic system at equi- 
librium depends on the distance from the rotation axis, an effect pointed out by Israel Q. However, the thermal 
equilibrium is characterized by a uniform temperature T as measured by the observer; another proof of this statement 
based on simple arguments can be found in Appendix B. 

The partition function Z^ defines an ensemble, that we can define as rotational grand- canonical where the angular 
velocity is fixed and the intrinsic angular momentum J can fluctuate. On the other hand, the grand-canonical 
partition function with fixed intrinsic angular momentum J defines an ensemble that can be called micro-rotational 
grand- canonical. The above definitions are inspired of the more familiar canonical- microcanonical duality, according 
to whether T or E are fixed. 

Pursuing the analogy, it is expected that these two ensembles become equivalent in the thermodynamic limit as 
far as the calculation of first-order moments of statistical distributions is concerned. In fact, neglecting the constant 
small factors multiplying Zj, it can be seen from ([5]) that logZ^ is essentially the Legendre transform of log Zj with 
respect to the total angular momentum J: 



log Z^ = log Zj 



logZj 



J-^log^ 



(14) 



Note that thermodynamic equilibrium is possible in a rotating cell because its acceleration is stationary 
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The saddle-point equation ([5]) for <p = iuj/T is simply the inverse Legendre transformation of (fl4 
j = Y-±L- j d 3 x / d 3 ptrO^(R,( W /T))(xxp)e^/V-( xx P)/ T 

t v^^y Jv J 



E 



Xj I d 3 x / d 3 p 



(27T) 



V 



'' trD^R^iw/r)) 



duj/T 



^ e -^/' r e"- (x><p) / T 



= L(w/T) + S(w/T) = log Z w (15) 

Therefore, the saddle point equation (|15[) relates angular velocity u> and angular momentum J and applies to relativistic 
as well as non-relativistic gases. It expresses the conservation of angular momentum, in the sense that the total angular 
momentum in the rotational grand-canonical ensemble on the right hand side, where the two contributions of orbital 
and spin can be clearly identified, equals the initial conserved total angular momentum J on the left hand side. It is 
worth pointing out that for macroscopic systems, the contribution of the orbital angular momentum is predominant. 
This can be seen from eq. (|15l) noting that the integral defining the total orbital angular momentum is scaled by a 
factor || xx p| with respect to the spin angular momentum. Since this is generally much larger than 1, for macroscopic 
distances and not too low temperatures, we usually have ||L|| ||S|| and the spin angular momentum can be neglected. 
Yet, for not too large systems, the spin angular momentum contribution can be sizeable. There is also a remarkable 
difference between them concerning their dependence on the size of the sytem: the spin angular momentum is properly 
extensive, i.e. it increase linearly with the volume, while the orbital angular momentum is not because of the factor 
|| x x p|| . The lack of extensivity of rotating systems makes the thermodynamic limit essentially irrelevant and dictates 
a finite-size treatment. From a completely equivalent point of view of the comoving frame, we can state that the 
rotating system is non-extensive because of the presence of a long-range centrifugal potential. 

Finally, it can be shown that in the rotational grand-canonical ensemble each microstate with total energy E, total 
charge Q and total angular momentum J has a probability (in the frame where (P) = 0): 

p(T, n, w) oc exp[-£;/T + [iQjT + u ■ 3/T] (16) 

and, consequently, rotational grand-canonical partition function can be written as: 

Z^ = tr{exp[{-H + fiQ + uj-J)/T}P v } (17) 

where J is the total angular momentum operator. This statement holds provided that ui/T C 1 so that terms of the 
order (u/T) 2 and higher can be neglected; only in this case is the formula (|17p equivalent to [[?]). Indeed, working out 
the trace (jTTJ) for an ideal gas involves the calculation of matrix elements such as: 

53<p, a\ exp[u; • J/T]P v \p, a) = ^<p, a\R*(iu>/T)P v \p, a) (18) 

a a 

where \p, a) are single particle states with four-momentum p and spin projection a (see Sect. 5 for extended discussion) 
and the rotation of an imaginary angle iuj/T has been introduced. The above single-particle trace has been calculated 
in ref. Q: 

Y,(pMR*{iu/T]Pv\p,°) = ^3 ^d 3 xe-(P- R -(-/ T )" 1 P)trZ? s (R li ^/T))(0|P^|0) (19) 

The vacuum expectation value of Py is immaterial and also becomes 1 in the large volume limit. Now, since: 

Rii(VO v = vcos?/> + (n x v) sin?/; + (1 — cos-0)v • nn (20) 

we have: 

x • (p R^iu/T)- 1 ?)] ~ i^x ■ (w x p) 

only if lo/T is sufficiently smaller than 1, so that cosine terms in (|20p can be approximated with 1 and sine terms 
with their argument. In this case, eq. (|19[) turns into, for large V: 

a\R*(iu,/T)P v \p, a) ~ ^3 J y d 3 x e^^' T tvD s '(Ra(iw/T)) (21) 



It is then straightforward to show that the final expression of (Tr?|) is just 0. 
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III. MOMENTS OF INERTIA 



In this section, wc will study more in detail the saddle-point equation (|5]) or its equivalent form (|15p . This is a 
vector equation in uj, whose solution depends, besides T and J, on the shape and size of the system. 

If the total angular momentum J is a simmetry axis of V, then J and uj must be parallel (see Fig. 1). This 
can be proved by applying an arbitrary rotation of an angle ip around J, let Rj (rp), to both sides of eq. (|15p . 
Since R(x x p) = Rx x Rp, defining new integration variables x' = Rx and p' = Rp, taking into account that 
uj ■ R _1 (x' x p') = Ruj ■ (x' x p') and being both momentum and spacial domains invariant by rotation around J, we 
can conclude that, if uj is a solution of the eq. Rj(rp)u> is also a solution for any rp. Therefore, Ri(rp)u> should 
coincide with uj for the equation to be well behaved. If this is the case, uj is parallel to J and the orbital term L(w/T) 
should be also parallel to uj. In this case, the eq. (|15|) essentially reduces to a scalar one: 

(2tt) 3 

j 



J = ^_^_ J d 3 x j d 3 ptr^(R,HT))J.(xxp)e-^V i -( xx ' > )/ T 

-ej/T e wJ-(xX P )/T fa) 



(2ir) s 



:j 



trD s i(Rt,(iu/T)) 



duj/T 



which can be solved numerically. 

Let us now introduce a length R related to the size of the system; it can be, e.g. the maximal distance of a point of 
the set V from the rotation axis. An analytic solution of the general equation (|15p can be obtained for small values 
of LoRp/T. In fact, under this circumstance, the exponential can be approximated as: 

exp[u; • (x x p)/T] ~ 1 + u) ■ (x x p)/T (23) 

This condition leads to different requirements for a non-relativistic and an ultrarelativistic gas (being p — 0(V mT) 
with m 3> T and p = 0{T) respectively), yet both imply that uiR <C 1. If approximation |23|) applies, the two terms 
L and S reduce to, at the first order in uj: 

L - E (^X d3x / d3p + ^ x x p)? • ( x x p ^ £j/T 

S - E (2^3 J v d3x J d3 P 3 6 3 T (24) 

By expanding the scalar and vector products in the integral of the orbital angular momentum and carrying out simple 
algebraic calculations, one obtains: 



3 (p 2 ) r 

u = E^)^ / d3x - x ^ Uk ( 25 ) 



where r 2 = X)i=i x f an d: 



^g^/d 3 pe"^ (26) 



is the mean density of the species j in the Boltzmann approximation and without angular momentum constraint, 
while: 

w=<h) X - ! W 1 l A ^'' IT (27) 

is the mean squared momentum of the particle j under the same approximations. 

The relation between L and uj (|25|) resembles the classical linear relation between angular velocity and angular 
momentum defining the inertia tensor I. Indeed, in the non-relativistic limit, this is precisely what one gets from (|25|) 
because (p 2 ) = ZmfT and, consequently: 

3 . 3 

Li = ^ E m o( n o) / ( r2 ^fc _ x i x k) w fc = E IikUJk ( 28 ) 

fc=l j V k=l 
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Hence the expression in (|25|) : 

= / d'x (r 2 <k- - x,x fc ) (29) 

turns out to be the right generalization of the inertia tensor for a relativistic gas. Like in the classical case, if the 
system (which is homogeneous for small lo) is symmetric around its rotation axis, then L is parallel to uj and the 
proportionality constant is just the moment of inertia with respect to that axis. 

Unlike in the classical case, for non-macroscopic systems there is also a potentially non-negligible spin contribution. 
In this case, the full relation between J and us reads: 




Z 



FIG. 1: For an axially symmetric system around the angular momentum direction, uj is parallel to J. Also shown the coordinates 
chosen to describe spectra in Sect. II Vl 



IV. SPECTRA AND ANISOTROPIES 



The impact of a hnite angular momentum on particle spectra in a relativistic gas was considered by Hagedorn Q 
who pointed out that a large value of the intrinsic angular momentum involves a peculiar anisotropy. Here, we will 
provide the most general expressions of the momentum spectrum of particles in an equilibrated rotating system with 
fixed (classical) angular momentum J. This can be obtained directly, for the jth particle species, from the partition 
function ([3]): 

^exp[(-H + fiQ)/T}PjP v 



/dnA_ J_ 



\d 3 p/ 



tr 



d 3 p 



where drij/d 3 p is the momentum spectrum operator. The above expression can be rewritten as the functional derivative 
with respect to a(p) of: 



z A a (p)] = tr \ ex P 



(-H + nQ)/T- 



|d 3 p«(p)^ 



d 3 p 



PjPv 



(31) 
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in a(p) = 0. Since, for an ideal gas: 



d 3 p 



the functional Zj[a(p)] can be simply obtained from eq. ((3J replacing £j/T with Ej/T+a(p). The functional derivative 
then turns out to be: 

(P) = A-^o g Zj[a( P )} = i- ^V \ 3 f d 3 xe^^ / d 3 0e^( J —P)tr^(R^(0)) 

\d J p/ 5a(p) «( P )=o Zj 8tH (27r) J „/ y J\d>\<ir 9 



x exp 



/ d 3 x /d 3 pe- Sfc / T e-^ ( xX P)trD s *(R^(^)) 



(32) 



Similarly to eq. ([3]), it is possible to make a saddle-point expansion of the above integral. If the orbital angular 
momentum of the particle is much lower than the total angular momentum, i.e. ||x x p|| <c J, the former can be 
neglected in the saddle-point equation which is then essentially the same as that for the partition function, eq. ([4]). 
This condition amounts to take exp[— icf>- (x x p)] as a constant factor and it is usually met except for momenta larger 
than J/R, R being the size of the system. In this case, the formula (f3"2"j) becomes, at the leading order: 



e -e S /T e u,-(x.xp)/T 



(33) 



which is, of course, the spectrum in the rotational grand-canonical ensemble in the Boltzmann limit. 

It is of some interest to derive the shape of the spectra in cylindrical momentum coordinates px = \Jpt + Vyi 
azimuthal angle if and rapidity y — 1/2 log[(p z + e)/(e — p z )\ when the angular momentum of the system is directed 
along the y axis, i.e. J = Jj and the system is rotationally symmetric around J. Because of this symmetry one has 
lj = uj] and the spectrum (|33f reads: 



, d 3 P / (2?r) 
In the new coordinates pT,f,y the spectrum reads 



^trZ^(Ra(*u;/T)) ^d 3 x- 



-ej/T w(xx.p)/T 



(34) 



\ prdprdipdy 



X j tiD s i(Ru,(iLj/T)) 
(2tt) 3 



d 3 x itlt cosh y exp [ — uix coshy/T + px • x x)±/T + rriT(^ x x)|| sinhy/T] 

(35) 

where wit = \J Pr + m j y || labels the longitudinal projection along the z axis and j_ the transverse projection onto 

the xy plane. 

With successive integrations of eq. (|3"5|) one obtains: 



/ dn j 
\pTdpxdip 



47T 3 



d 3 x 



1 



1 — (U3 X x) 



and: 



= totKi (rriT\J^ — (u; x x)^/Tj exp[py • (a; x x)j_/T] 

(36) 



/ dnj \_ XjtrD s ^R^(iuj/T)) 
\p T dp T / 2n 2 



d 3 x 



v Jl-(ux x)j 



: m T K! (mryjl - (w x x)|/r) I (pt|I(" x x) ± ||/T) (37) 



where Ki and Io are McDonald and modified Bessel functions respectively. It is also possible to calculate the mean 
multiplicity by integrating (|37j) . which yields: 



^taD s '(M»w/T)) / d 3 X i f T K 2 (mVl - ||u> x x|| 2 /T 



(38) 



The formulae (f3"5]) - (f3"8")) yield as limiting cases for <jj = the well known ones for the relativistic ideal Boltzmann gas. 
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It is worth pointing out that the spectrum ((35)) features a non-trivial dependence on the momentum azimuthal 
angle ip strongly sensitive to the ratio lo/T. Introducing cylindrical coordinates r, $,z for the vector x and taking 
into account that u) — u>j, eq. (|36p can be written as: 

/ drii \ A,trD^ (R^(iw/T)) f , q 1 / , -, , \ 

( 2 -a ) = s / d 3 x =m r K! (rriTv/l- wrcos$ 2 /T) exp[p T coz cos ^/T 

XPTdpTd^/ 47r J J y ^Jl - ( W rcos$) 2 V / 

which shows that particles emitted along the momentum axis (with ip = tt/2) have, on average, a lower momentum 
than those emitted orthogonally to it. Classically, this can be understood as an effect of centrifugal force. 

The (|39[) is a periodic function in the momentum azimuthal angle tp, thus it can be Fourier expanded. Borrowing 
a notation commonly used in relativistic heavy ion physics, we will denote by 2vk the coefficients of this expansion, 
i.e. setting: 

/ \ ^ / \ 2 ' 2k (40) 

\p T dp T dip J 2tt \p T dp T I Z-j 2 



with vq = 1 we obtain: 



JO \pTdpT<itp/ r 

= T~~ v (41) 



Jo 



and, by using |39 



Jv d3x ^f TO ^ K i f^T^ 1 - W^W/t) hkiPTZtu/T) 

w 2fe = 7 — x (42) 

J v d 3 x ^_ ( j co9$)a "irKi (mry/l- (cur cos®) 2 /T) I (ptzuj/T) 

which shows that the Fourier coefficients are bound between and 1. For low px, the ratio between modified Bessel 
functions is such that l2fc/Io oc p?^, while at large px values, the ratio of modified Bessel functions tends to 1 and, as 
a consequence, all V2k — * 1. 



V. POLARIZATION 



If a gas has a net electric charge, an imbalance in the mean multiplicities of positive and negative particles is 
implied. Similarly, if a gas at thermodynamical equilibrium has a non-vanishing angular momentum (hence it is 
rigidly rotating, as we have seen in Sect. [TTJ) particles should have a net polarization along the direction of angular 
momentum. From the viewpoint of the comoving observer in the rotating frame, the hamiltonian in the rotating 
frame has a spin-rotation coupling term 0, [13] 2 - For the non-relativistic case, the calculation is straightforward as 
the polarization vector is the same in the observer and in the comoving frame, hence it can be carried out in the latter 
where the single-particle hamiltonian reads (lpj : 

h xo t = h ohs - uj ■ j (43) 

being j the total (orbital + spin) angular momentum operator of the particle. In fact, the problem we want to solve 
here is to calculate the polarization vector for a relativistic gas. 

This problem is more difficult than for charge, energy or momentum because angular momentum is not a generator 
of an abelian group and this adds some complication. Particularly, from the point of view of statistical mechanics, 
one has to deal, as we will shortly see, with a spin density matrix with non- vanishing off-diagonal elements, unlike for 
momenta. From eq. (|17p. the density operator in the rotational grand-canonical ensemble reads: 

% = ^cxp[(-i/ + /iQ + o;- J)/T]P y (44) 



2 This phenomenon is well known for radio waves [3 
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while in the micro-rotational grand-canonical ensemble, from eq. ([3]) : 

P. J = T^exp[(-iT + /4)/T]PjP v (45) 
A/ 

Pj being the quantum projector onto states with definite angular momentum. As argued in Sect. [Til the two ensembles 
are equivalent for large J and V provided that oj/T is sufficiently smaller than 1. Henceforth, we will confine our 
attention to this case and we will work in the rotational grand-canonical ensemble described by (|44|) . 

In the Boltzmann limit of the ideal relativistic quantum gas, all particles can be handled as independent distin- 
guishable objects. Thus, the density operator (|4"4"]) factorize and we can calculate the polarization by considering the 
single-particle density operator: 

p w = — cxp[{-h + pq + u • j)/T]P v (46) 

where 2 W is the single particle partition function: 

z u = |^d 3 x J d 3 p e-z/ T e"<*^ T trD s (R^Lo/T)) (47) 

ensuring the correct normalization trp w = 1 (this in fact requires ui/T <C 1, see discussion at the end of Sect. UT|) . 
The density operator restricted to the spin degrees of freedom is a function of momentum and reads: 

MP) = I ^ F/y d3xe-^-(xx P )/Ttr Jg (R,^/T)) + ^ + " ^ 

The polarization is, by definition, the trace of the density matrix p(p) multiplied by a suitable spin operator. In 
non-relativistic quantum mechanics, this is obviously the spin vector operator S and the polarization can be obtained 
straightforwardly: 



II = tr[SA»] = ^" S -~ S '~~ — C (49) 



which turns out to be momentum- independent. In relativistic quantum mechanics the proper generalization of the 
spin angular momentum is the Pauli-Lubanski vector multiplied by l/m: 

% = -(1/2) ^v P J vp P° (50) 

which in fact reduces to the spin operator in the particle rest frame; J vp are the generators of the Lorentz group. The 
polarization is thus a four-vector II: 

n(p) = -tx p (Wp(p)) (51) 

m 

where the trace is to be calculated by summing only over spin degrees of freedom keeping the four-momentum p fixed. 
This four-vector Tl(p) has vanishing time component in the particle rest frame, as it is apparent from (|50[) . Before 
working out eq. (|5ip . we should introduce some important notions about the construction of physical states. We will 
stick to the notation of ref. [ll| . 

As it is well known, the Pauli-Lubanski vector fulfills the commutation relations: 

[%,P V }=0 (52) 

pa 

W P = 

Hence, if the ket \p) is an eigenvector of P, so is The restriction of W to the eigenspace labelled by four- 

momentum P is defined as W(p). Since W(j>) ■ p = 0, this four-vector operator can be decomposed onto three 
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orthonormal spacelike four- vectors ni(p), n^ip), n%(jp) which form a basis of the Minkowski space with the unit vector 
P = p/\fp 2 - 

3 

W{p) = J2w i {p)n i (p) (53) 

i=l 

It can be shown that the operators: 

Si(p) = Wi(p)/m (54) 

form an SU(2) algebra and are the actual relativistic generalization of the spin angular momentum. The third 
component S^p) can be diagonalized along with S 2 — —W 2 /m 2 which is a Casimir of the full group 10(1, 3)^, with 
corresponding eigenvalues a and S(S +1). With a suitable choice of rii(p), i.e.: 

niip) = ]?]<* [p] = R3(p)R 2 (fl)L 3 (0 (55) 

ej being the unit vectors of spacial axes and [p] being a Lorentz transformation bringing the timelike vector po = 
(m, 0, 0, 0) into the four-momentum p with polar coordinates £, 9, tp; the eigenvalue A has the physical meaning of the 
component of intrinsic angular momentum in the rest frame along the direction of particle momentum p. Thus, with 
the choice (|5"5j) . A is the helicity in the rest frame and J is, by definition, the spin of the particle. Since, from eqs. (|54p . 
([53")) and (|50[) Si(po) — Ji, the spin operators in the rest frame coincide with the generators of the rotation groups. 
Finally, single particle states will be written as \p, a) with: 

P\p, a) = p\p, a) and S 3 (p)\p,a) = a\p,cr) (56) 

and normalization: 

(p,a\q,T) =5 3 (p-q)5 aT . (57) 
while the transformation of a state \p, a) under a general Lorentz transformation A reads: 



A\p,a) =^|Ap,r)^ CT ([Ap]- 1 A[p])y^ . (58) 

We are now in a position to develop eq. (JHH). By using (fS"3")) and we get: 

1 _ 3 

H W = -T,^ a \ W \P' a '^ a '\MP)\P^) = ^^<P,v\Si(p)\py)<py\fa(j>)\p,<r)ni(p) (59) 

Two different matrices appear in the above equation. By definition: 

{p,a\S i (p)\p,a , )=D s arT ,{J i ) (60) 

are the matrices of the SU(2) generators Ji in the representation labeled by the particle spin S. Furthermore, according 
to 



{p,o'\Mp)M = ^ j v d3x e ., xxp ) /T tlD s {R ^ lu}/T)) <p. °'\ e ^ -VVPvM (6i) 

As to the rightmost matrix element, we will restore the vector <p = ilj/T and calculate it for real <p (or imaginary u>), 
then making an analytic continuation to imaginary cp, the same way we did at the end of Sect. [TTJ Therefore: 

(p, CT '|exp[u;-j/T]P y |p,a) = (p, a'\ fflq>[-*0 • j[P v \p, a) = (p, tr'|R^)Py |p, a) (62) 

The latter matrix element has been calculated in ref. [|[ 3 : 

(p )( 7'|%^)P v b,(7) = |d 3 p' <S 3 (R^) P '- P )^ |d 3 xe ix CP-P') 

1 



x-(^(b]- 1 R^(0)[p]) + ^(b]tR^(0)[p]t-i ) )^(O|P y |O) 



(2tt) 3 



jf d 3 xe*-fr- R *« _lp >±(l? s (^ (63) 



3 It can be easily obtained from eq. (55) in the reference 
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In this equation, [p] is now meant as an element of the universal covering group of Lorentz group, i.e. SL(2,C), so that 
the notation [p]^ becomes meaningful as well as its finite-dimensional representation matrix D s ([p\) corresponding to 
the Lorentz group representation usually labelled as (5, 0) (the (0,5) being the D s ([p])^ -1 ). The factor (0|Py|0) is 
immaterial and also becomes 1 in the large volume limit as Py — > I. 

If 4> <C 1, as it was supposed to be at the beginning of this section, then: 

ax- (P - R$,(0) _1 p) - • (x x p) 

and we can rewrite eq. (|6"5|) as: 

{p,a'\R^)P v \p,a)^^ ^d 3 xe-^( xx P)i(^(H- 1 R^(0)b])+ J D s (b]tR^(0)[ P ]t-i))^ CT (64) 
We can now make an analytical continuation to imaginary <f> or real lo and obtain, by using (frITj) and (|M)) : 

<py\p u (p)\p,<r) = 2trD s {R \ {ioj/T)) (^ g (M" 1 ^(^/r)M) + ^ 5 (M t Ra(Wr)M t " 1 ))^ (65) 

This matrix is hermitian and has trace equal to 1, as required for a good density matrix. Hermiticity can be shown 
by taking advantage of a remarkable feature of SL(2,C) representation [i~l| : 

D s (A' f ) = D s (Ay 

and taking into account that D (R&^w/T)) is hermitian. 

Plugging (|65|) into (|59|) . using ([60]) and carrying out the momentum integration, the final expression of the polar- 
ization four-vector for a particle of momentum p in the observer frame is finally achieved: 

n{pT = 2tr^( R l( w /r)) E tr [^WW^WWI+^WW^MIP] 1 " 1 )] «i(P) 

3 

J> [^(JiW-^fiw/TJH) + i? s (J,:H t R^^/T)[p]t- 1 )] ([p] ei f (66) 



2tr^ s (R^(^/T)) ^ 

The proper polarization four-vector, in the particle rest-frame, has components ([p] -1 !!^))^, that is: 

MpT = 2trgS(R ^ (Wr)) E tr ^ s Wb]" 1 R^(^/r)b]) + £> s (J < b] t ^(^/r)[p] t - 1 )] sf (67) 

which has vanishing time component, as required. 

The equations (f6"6")) and (j6"7| are the general analytical expression of the polarization for a particle with spin 5. 
We will now develop them for the most interesting cases of spin 1/2 and spin 1. Before doing this, we first observe 
that the Lorentz transformation [p] to calculate the polarization can be chosen arbitrarily. In fact, if [p]' ^ [p] also 
transforms the timelike vector po = (to, 0,0,0) into p, then [p] _1 [p]' is a pure rotation R as it leaves po invariant. 
Hence, the polarization four-vector II(p) defined with [p]' becomes: 



u( P ) cx ^^[^(j^r^^^/^Mo+^^JiW'tR^^/^M't- 1 )] (\py ei ) 
»=i 

3 i 

= ^ -tr [D^JiR-^-^iiuj/TMR) + D s (J.R-^R^/T) M t_1 R)] (MRe<) 

i=l 
3 i 

= Eo tr [^ S (R^R" 1 W" 1 R^(^/r)[p]) + ^ s (RJ l R- 1 [p] t R li (Wr)b] t - 1 )] (tp]Re0 



3 



£ J tr [D s (J Ret [p] - 1 R^ (iw/T) [p] ) + # s ( J Res [p] t R* (tw/T) [p] t~ 1 )] ( [p] Re,) 



2 

i=l 



(68) 



where we used the unitarity of R. The last expression in (|68p is apparently equal to the one in (|66[) . 
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A. Spin S=l/2 

Choosing for [p] a pure Lorentz boost, we have fill ]: 

^ 2 (M) = ^ ( |p])t = m + e + g -P 

y 2m{m + ej 



and 



^ (b]) -l =jD V2 ([p]) t-l = 



m + e — cr • p 

y / 2m(m + e) 



where <r are the Pauli matrices. Therefore: 

tr [^(J^l^Miw/TM 
whereas, because of (|69|) and ((70)) : 

tr [D^iJ^Roiiw/T)^- 1 ) 
i.e. it is obtained from (|7"Tj) by reflecting p. Since 



= tr 



cr,: to + £ - er • p /rTn x"i + e + cr-p 

Rc{iuj/T) — ) 



2 y/2m(m + e) 



y/2m{m + e) 



tr 



CTim + e + cr-p m + e - <r • p 



2 y/2m(m + e) 



^J2m(m + e) 



(69) 



(70) 



(71) 



(72) 



D ' [Rmiuj T)) — I cos icr ■ u) sin — = I cosh h cr ■ u> smh — 



we have, for (fTTj) 



m + e — cr • p 

y/2m(m + e) 



= 1 cosh — — + smh — — 
2T 2T 



1 cosh h cr ■ uj smh — — — 

2T 2T/ 1 v /2m(m + e) 



to + £ i . „ . UJ ■ per ■ p — cr ■ (p X (ujx p)) 
— cr • u> H cr ■ wxp — ; r 

2m m 2m(e + to) 



(73) 



(74) 



Now we have to take the sum of the two traces (f7"Tj) and ([72]) implying that all terms which change sign in a reflection 
of the momentum p in (|74p vanish. Therefore, we are left with: 



2tr 



T 1 W • 1 U 

1 cosh \- sinn — 

2T 2T 



m + £ „ o> ■ per -p-cr ■ (p x (w x p)) 
cr ■ uj 



2m 



2m(e + m) 



sinh • 



2T 



to + e „ w • ppi - (p x (d> x p))j 

Wi 



m(e + m) 



Taking into account that: 



trD^ 2 (R^(iuj/T)) = 2 cosh — 



(75) 



(76) 



we can obtain the polarization vector in the rest-frame putting (|75|) and (I76p into eq. (|6T[) : 



^ 1 , ^ 
Lin = — tanh — 

4 2T 

1 , w 
= — tanh — 
4 2T 

1 , 

= — tanh — 

2 2T 



to + e „ d> • pp — p x (c2> x p) 
a; - 

to m(s + to) 

to + e „ a> • pp — p 2 a> + Cj ■ pp 

UJ 



m(e + to) 



uj pp 



777 m(s + to) 



£ 

UJ 



(77) 



Therefore, unlike in the non-relativistic case, the polarization vector has a component along particle momentum in the 
observer frame. This effect is owing to the vector nature of the polarization; indeed, the components of polarization 
in the observer frame ((66]) can be obtained with a general Lorentz boost [13]: 



n = n n 



7 + 1 



/3-n 0/ 3 



n° = 7/3 • n 



(78) 
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with (3 = p/e and 7 = e/m. Therefore, by using ([77)1 , eq. (|78|) becomes: 



1 , w e „ 

II = - tanh — — — u> 
2 2Tm 



n° = - tanh - 

2 2T m 



(79) 



Hence, the polarization three- vector is aligned with the angular velocity in the observer frame and not in the particle 
frame. Yet, there is also a time component which vanishes in the non-relativistic limit p/m — > 0, where the (|79|) 
correctly yields (|19")) . The remarkable difference with respect to the non-relativistic case (pit?)) is that polarization now 
depends on momentum. The longitudinal component, along p, is: 



n • p = - tanh — p • w 



while the component along the rotation axis turns out to be: 



n ■ u> 



1 , w 

— tanh — 

2 2T 



£ _ P 2 (^ • P) 2 
m m(e + to) 



(80) 



(81) 



which shows a very interesting feature: the polarization is maximal for particles with momentum orthogonal to 
the rotation axis. Furthermore, the polarization increases with energy, being proportional to e/m. Of course this 
behaviour cannot go on indefinitely because polarization cannot exceed 1/2 in any direction; since w/T -C 1, it is 
seen from eq. (|81|l that something must happen when e ~ 2mT/u>. Indeed, we have pointed out in Sect. IIVI that 
when momentum is of the order of J/R, the saddle point expansion of partition function at fixed J can no longer be 
independent of p, and uj becomes in fact a function of p, that is to(p), for the spectrum (l32|) expansion. This has some 
impact on the polarization vector: at some large momentum the dependence of lo on p should restore the natural 1/2 
bound. 



B. Spin S=l 



The calculation of the polarization for massive spin 1 particles is more involved than for spin 1/2 and we have carried 
it out by performing explicitely the multiplication of matrices of SL(2,C) D 1 representation. Since eventually one has 
to calculate traces, the choice of the basis for the matrices is arbitrary and we have written [p] = ^3(^)^2(^)^-3(0 ' 
being cosh£ = e/m, in the cartesian basis for the D 1 representation space: 



D x (R 3 (tp)) 



cos <p — sin <p 
sin ip cos (p 
1 



D\R 2 (9)) 




D\L 3 m = expH^ 3 ] = exp[£J 3 ] = D^R^)) 



D x (R C) (iu;/T)) = exp[u;/r • D X (J)] 



i)jk 



cosh£ — isinh^ 
i sinh £ cosh £ 
1 



Multiplying the matrices in eq. (fB"2")) and calculating the traces, according to eq. (|S7)). one gets: 



2sinh(w/T) 
2cosh(w/T) + 1 



w ■ pp 



m m(e + m) 



(82) 



(83) 



where the denominator 2 cosh(o;/r) + 1 is the trace of the matrix D 1 (R l ~,(ioj/T)). Therefore, the kinematical structure 
of the polarization vector is the same as in the spin 1/2 case, which is a reasonable outcome; in fact, the properties 
of Lorentz transformation of a polarization vector should not depend on the particle spin itself. 

In the case of spin 1 it is also interesting to calculate the fraction of the 00 component of the density matrix (|65l) . 
The calculation can be done quickly by noting that: 

Puooif) = trP 3 p = J2(P^\ p 3PUp)\P,o) (84) 

a 

P3 being the projector onto the state \p, 0). Written in the cartesian basis the matrix corresponding to P3 is simply: 

(85) 
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because of the choice of z axis as quantization axis, according to (|56| Therefore, by plugging eq. (|65|) into eq. I|84p: 

1 



PwOo(p) = 



2trL> 1 (R^(«w/T)) 



tr [PsDWpr^iw/TM) + P.DWp^R^/T)^- 1 )} 



yielding: 



Puj oo (p) 



1 



2cosh(w/T) + 1 



For small w/T, we have: 



Puj oo (p) 



cosh^/T) + % ; (1 - cosh(w/T)) 



p 2 w 2 



1 l-3(p-a>) 2 u 2 
3 + 18 T 2 



(86) 



(87) 



Discussion 



In both spin 1/2 and spin 1 cases, we have seen that the proper polarization vector has the same structure, so that 
we can easily generalize the above results to any spin: 



, ne 



j/T 



V 



e A u> ■ pp 

— IjJ 

m m(e + to) 



(89) 



This non- vanishing polarization means that spin states are not evenly populated in an equilibrated thermodynamical 
system with finite, macroscopic, angular momentum. Of course reaching equilibrium for the spin degrees of freedom 
implies that a small interaction should exists involving particle spin, e.g. a spin-orbit coupling. Yet, for the ideal 
relativistic gas, this interaction is assumed to be negligible in comparison with the pure kinematical effect of angular 
momentum conservation. This can be rephrased in the rotating frame as follows: the actual interaction involving 
spin is negligible with respect to the spin-rotation coupling. Since a macroscopic thermodynamical system can be 
conceptually divided into elementary fluid cells, and each cell participating in the rigid rotation is an accelerated 
system, for our result to be consistent with locality, one has to conclude that, in general, acceleration involves a 
polarization expressed by (|89[) where u) is to be interpreted as a local vector field involving local quantities such as 
velocity and acceleration of the cell. 
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APPENDIX A - ENTROPY FOR A SPINNING SYSTEM 



The entropy of a system with finite angular momentum J can be calculated easily from the expression of the 
probability of a state with energy E, vanishing momentum, charge Q and fixed J in the grand-canical ensemble: 

p=± exp [-E/T + »Q/T] (90) 

A7 



so that, by using (|14p 

stateswithfixedJ 

which is the known expression of entropy for a rotating system. The logarithm of the partition function can then be 
identified with the integral of the pressure: 

logZ w = l J d 3 xp(x) (92) 
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where, unlike in familiar cases, the pressure is not uniform due to rotation. From the entropy expression: 

TS = U - fx(Q) - u ■ J + / d 3 x p(x) (93) 

Jv 

we can derive the relation Q 



dS_ 
83 



= (94) 

V,U,Q T v 1 



APPENDIX B - EQUILIBRIUM CONFIGURATION OF A RELATIVISTIC ROTATING SYSTEM 

We will now show that a relativistic macroscopic system with finite angular momentum at equilibrium must be 
rigidly rotating by generalizing an argument by Landau [5J for non-relativistic system. 

Let us consider a generic isolated hdyrodynamical system with velocities Vj i being the label of hydrodynamical 
cells. In order to find the equilibrium configuration the entropy should be maximized with the constraint of energy, 
momentum and angular momentum conservation. Therefore, we have to find the extremum points of: 

J2S i -P(^2E i -E ) + /3-J2 p i + ^-(E / Xi x P « - J ) ( 95 ) 

i i i i 

where we have taken the total momentum vanishing, i.e. we are working in the system's rest frame; Eq is the total 
energy and J the total angular momentum; /3, (3 and lj/T are Lagrange multipliers enforcing the conservation laws. 
Since the entropy is a relativistic invariant, it can only depend on the mass of the cell, i.e. Si = Si(y/ Ef — Pf ). In 
order to find the equilibrium configuration, one has to maximize (|95p with respect to each Pj and E^. Thus: 

dSi Ei dSi 

dEr^dwr^ yi (96) 

and: 

Taking into account that dSi/dMi is the inverse of the proper temperature Ti of the ith cell by definition, eq. (|96|) 
implies that: 

— = constant = 8 = — 
Ti T 

where T is defined the global temperature of the system. Plugging (|9"8|) into (f9"T|) we get: 

7,-v,- V,- uj x Xj 

— = — =(3+ 99) 

For the total momentum to vanish, the vector (3 should be as well and we are left with: 

Vj = uj x Xj (100) 

that is a rigid rotation around the axis uj. To show that T = 1/8 is in fact the global temperature of the sytem, one 
has to consider that at equilibrium, using formula (|91[) in Appendix A for the entropy of the cell i: 

i i i i 

M being the mass of the system, because P = 0. Therefore: 

dS _ 1 
DM ~ T 

showing that T is the actual temperature of the system. The relation l|98p can be rewritten as, by using (|100p : 

T = 7Vl-||u;xx l || 2 (102) 
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in accordance with eq. (jlip . The involved physical meaning is as follows: the temperature measured by a thermometer 
at rest in the observer frame, is lower than that measured by a comoving thermometer, i.e. at rest in the cell frame. 
This is a pure relativistic effect which has no correspondance in classical thermodynamics. 

The above argument can be extended to include the intrinsic angular momentum contribution Ji of each cell to J. 
Eq. ([95]) now reads: 

Si(y/E?-Pl Ji) - f3(J2 Ei - E ) + (3 ■ V* + P<» ■ (J2 Xi x P ' + 3 > - J ) ( 103 ) 

i i i i 

The conclusions are similar, with the additional condition: 

dSi u> 
dJ~l = ~T 

in accordance with eq. (|94l) . 
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